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Abstract 



A random-matrix theory is presented which shows that breaking time-reversal 
symmetry by itself does not suppress the weak-localization correction to the 
conductance of a disordered metal wire attached to a superconductor. Sup- 
pression of weak localization requires applying a magnetic field as well as 
raising the voltage, to break both time-reversal symmetry and electron-hole 
degeneracy. A magnetic-field dependent contact resistance obscured this 
anomaly in previous numerical simulations. 
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Typeset using REVT^ 



Weak localization is a quantum correction of order e^/h to the classical conductance of 
a metal.SH The word "localization" refers to the negative sign of the correction,! while the 
adjective "weak" indicates its smallness. In a wire geometry the weak-localization correction 
takes on the universal valueB 5G = — |e^//i at zero temperature, independent of the wire 
length L or mean free path £.@ The classical (Drude) conductance Go — i^^/L) e'^/h is much 
greater than 6G in the metallic regime, where the number of scattering channels N ^ L/£. 
Theoretically, the weak-localization correction is the term of order in an expansion 
of the average conductance {G) = Gq + 6G + 0{N~^) in powers of N. Experimentally, 
SG is measured by application of a weak magnetic field B, which suppresses the weak- 
localization correction but leaves the classical conductance unaffected.^ The suppression 
occurs because weak localization requires time-reversal symmetry (T). In the absence of T, 
quantum corrections to Go are of order A^~^ and not of order N^. As a consequence, the 
magnetoconductance has a dip around B = of magnitude 6G and width of order Be (being 
the field at which one flux quantum penetrates the conductor). 

What happens to weak localization if the normal-metal wire is attached at one end to a 
superconductor? This problem has been the subject of active research.Lryi The term Go of 
order is unaffected by the presence of the superconductor.! The 0{N^) correction 6G is 
increased but remains universal ,00 



6G 



2 - 871-^) eyh ^ -1.19eyh. 



(1) 



In all previous analytical work zero magnetic field was assumed. It was surmised, either 
implicitly or explicitly^ that 5G = in the absence of T — but this was never actually 
calculated analytically. We have now succeeded in doing this calculation and would like to 
report the result, which was entirely unexpected. 

We find that a magnetic field by itself is not sufficient to suppress the weak-localization 
correction, but only reduces 6G by about a factor of two. To achieve 6G = requires in 
addition the application of a sufficiently large voltage V to break the degeneracy in energy 
between the electrons (at energy eV above the Fermi energy Ep) and the Andreev-refiected 
holes (at energy eV below Ep). The electron-hole degeneracy {T>) is effectively broken 
when eV exceeds the Thouless energy Ec = hvpi/L"^ (with vp the Fermi velocity). Weak 
localization coexists with a magnetic field as long as eV <^ Ec. Our analytical results are 
summarized in Table |. These results disagree with the conclusions drawn in Ref. |] on 
the basis of numerical simulations. We have found that the numerical data on the weak- 
localization effect was misinterpreted due to the presence of a magnetic-field dependent 
contact resistance, which was not understood at that time. 



-6G [eyh] 


T 


no T 


V 


2 - S/tt^ 


2/3 


no V 


4/3 






Table ||: Dependence of the weak-localization correction 6G of a normal-metal wire attached 
to a superconductor on the presence or absence of time-reversal symmetry (T) and electron- 
hole degeneracy {V). The entry in the upper left corner was computed in Refs. [10], |Tl. 
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Starting point of our calculation is the general relation between the differential conduc- 
tance G = dl/dV of the normal-metal-superconductor (NS) junction and the transmission 
and reflection matrices of the normal region,! 

G = (Ae^h) tTm{eV)m\eV), (2a) 
m(e) = t'{e) [1 - a{ey{-e)r{e)]-^ t*{-e), (2b) 

where a{e) = exp[— 2i arccos(£/A)]. Eq. (|^) holds for subgap voltages V < A/e, and 
requires also A ^ Ep (A being the excitation gap in S). We assume that the length L of 
the disordered normal region is much greater than the superconducting coherence length 
^ = (^fi?£/A)^/^. This implies that the Thouless energy <^ A. In the voltage range 
V ^ Ec/e we may therefore assume that eV <C A, hence a = —1. The N x N transmission 
and reflection matrices t, t', r, and r' form the scattering matrix S{e) of the disordered normal 
region (A^ being the number of propagating modes at the Fermi level, which corresponds to 
e = 0). It is convenient to use the polar decomposition 

Vt iVr 

Here vi, V2, wi, and W2 are N xN unitary matrices, T is a diagonal matrix with the N trans- 
mission eigenvalues Tj G [0, 1] on the diagonal, and R = 1 — T. Using this decomposition, 
and substituting a = — 1, Eq. (pEj) can be replaced by 






m{e) = ^T{e) \l + u{e)^ R{-e)u*{-e)^R(^] 



■n -1 



■u{e)^T{-e), u{e) = w^{e)wl{-e). (2c) 

We perform our calculations in the general framework of random-matrix theory. The 
only assumption about the distribution of the scattering matrix that we make, is that it 
is isotropic, i.e. that it depends only on the transmission eigenvalues.lll In the presence of 
T (for B <C i?c), S = S^, hence Wi = wj- (The superscript T denotes the transpose of a 
matrix.) If T is broken, wi and W2 are independent. In the presence of V (for eV <^ Ec), 
the difference between StfV) and S{—eV) may be neglected. If V is broken, S{eV) and 
S{—eV) are independent .li^ Of the four entries in Table |, the case that both T and V are 
present is the easiest, because then u = 1 and Eq. (^ simplifies tS 

G={Ae'/h)Y.nTl{2-T^)-\ (3) 

The conductance is of the form G = J2n f{Tn), known as a linear statistic on the transmission 
eigenvalues. General formulasSIll for the weak-localization correction to the average of a 
linear statistic lead directly to Eq. (|I]). The three other entries in Table |, where either T 
or V (or both) are broken, are more difficult because G is no longer a linear statistic. We 
consider these three cases in separate paragraphs. 

(1) V, no T — Because of the isotropy assumption, wi and W2, and hence u, are uni- 
formly distributed in the unitary group U{N). We may perform the average (■ ■ ■) over the 
ensemble of scattering matrices in two steps: (■ ■ ■) = ((■ ■ ■)u)t) where (■ ■ ■)„ and (■ ■ ■)t are, 
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respectively, the average over the unitary matrix u and over the transmission eigenvalues Tj. 
We compute (■ ■ ■)„ by an expansion in powers of N~'^ . To integrate the rational function 
(0) of u over U{N), we first expand it into a geometric series and then use the general rules 
for the integration of polynomials of u. The polynomials we need are 

p,q=0 

Mpg = {tTT{uVRu*VRyuTu\VRu^VRu^y)u. (4b) 
Neglecting terms of order A^~^, we find 



Nrfil-T^fP if p = g, 
n{rf + Ti - 2t^){1 - ri)P+'?-i - 2 min(p, q) 
X T-fi^f - ^2)(1 - ^1^+9-2 if |p - g| odd, 
else. 



where we have defined the moment = N ^ X)j . The summation over p and q leads to 

^ _ Ati-2tI + 2tI-At2 

4e^ 2 — Ti Ti[2 — Tiy 

It remains to average over the transmission eigenvalues. Since is a linear statistic, we 
know that its sample-to-sample fluctuations 5Tk = — (rfe) are an order smaller than 
the average .0 Hence 

{f{n))r = f{{n))[l + 0{N-% (6) 

which implies that we may replace the average of the rational function @ of the r^'s by the 
rational function of the average (rfe). This average has the 1/N expansion 

(rfc) = (rfe)o + 0(iV-2), (7) 

where (Tfc)o is 0{N^). There is no term of order A^~^ in the absence of T. From Eqs. (iD-(0) 
we obtain the expansion of the average conductance, 

h N{n)o 4(ri)o-2(n)g + 2(n)g-4(r2)o 

Eq. (P) is generally valid for any isotropic distribution of the scattering matrix. We apply 
it to the case of a disordered wire in the limit oo, £/L — > at constant Ni/L. The 

moments (Tfc)o are given byll 

(ri)o = i/L, (r2)o = f^/L. (9) 

Substitution into Eq. (||) yields the weak-localization correction 5G = —^e^/h, cf. Table |. 

(2) T, no P — In this case one has u'^{—eV) = u{eV) and u{eV) is uniformly distributed 
in U{N). A calculation similar to that in the previous paragraph yields for the average over 
u: 
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^ : {G}u = iVri+r- 

2ri+r2_ - T^^rt - r^+rf_ - r^+rt - rl^r^_ + t^^tI_ + tI_^t^_ , (10) 



X 



where we have abbreviated = Tk{±eV). The next step is the average over the transmis- 
sion eigenvalues. We may still use Eq. (||), and we note that {Tk{s)) = (r^) is independent of 
e. (The energy scale for variations in {Tk{e)) is Ep, which is much greater than the energy 
scale of interest Ec.) Instead of Eq. (^ we now have the 1/A^ expansion 

{Tk) = {Tk)o + N~'6n + 0{N-^), (11) 

which contains also a term of order A^~^ because of the presence of T. The 1/A^ expansion 
of {G) becomes 

h N{t,)o , 26n , 2(ri)g-2(ri)|;-2(r2)o + 2(ri)o(r2)o , 

- + (2 - (r,)o)^ + (r,)o(2 - (n),)^ + ^^"^ ^^'^ 

For the application to a disordered wire we use again Eq. (^ for the moments (Tfc)o, which 
do not depend on whether T is broken or not. We also need 6ti, which in the presence of 
r is given byi 5ri = -i. Substitution into Eq. (|12D yields 6G = -^e^/h, cf. Table |. 

(3) no T, no P — Now u{eV) and u{—eV) are independent, each with a uniform distri- 
bution in U{N). Carrying out the two averages over U{N) we find 

A (G). = . (13) 

The average over the transmission eigenvalues becomes 

where we have used that 5ti = because of the absence of T. We conclude that 6G = in 
this indicated in Table DO. 

This completes the calculation of the weak-localization correction to the average con- 
ductance. Our results, summarized in Table |, imply a universal B and ^-dependence of 
the conductance of an NS microbridge. Raising first B and then V leads to two subsequent 
increases of the conductance, while raising first V and then B leads first to a decrease and 
then to an increase. 

So far we have only considered the 0(A^°) correction 6G to (G) = Go + 6G. What about 
the 0{N) term Go? From Eqs. (D, (pj), and (|l|) we see that if either T oi V (or both) 
are broken, 



G<,= '-^.!^ = (2eVh)N{l + L/(y\ (15) 



;2 

h 2- (ri)o 

In the second equality we substituted^ {ti)o = (1 -|- L/£)^^, which in the limit ijL — > 
reduces to Eq. (01). If both T and V are unbroken, then we have instead the resultli^ 
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Go = (2eV/i) N[l + L/i + 0{i/L)]-^ 



(16) 



The difference between Eqs. ([15|) and (0) is a contact resistance, which equals h/ANe^ in 
Eq. (p!5| ) but is twice as large in Eq. (|1^). In contrast, in a normal- metal wire the contact 
resistance is h/2Ne^, independent of B or V . The B and V^-dependent contact resistance 
in an NS junction is superimposed on the B and V^-dependent weak-localization correction. 
Since the contribution to (G) from the contact resistance is of order (e^//i)A^(£/L)^, while the 
weak-localization correction is of order e^//i, the former can only be ignored if N{i/Ly <C 1. 
This is an effective restriction to the diffusive metallic regime, where i/L <^ 1 and Ni/ L ^ 1. 
To measure the weak-localization effect without contamination from the contact resistance if 
N{i/LY is not <^ 1, one has two options: (1) Measure the 5-dependence at fixed V ^ Ec/e; 
(2) Measure the ^-dependence at fixed B ^ B^. In both cases we predict an increase of the 
conductance, by an amount |e^//i and |e^//i, respectively. In contrast, in the normal state 
weak localization leads to a 5-dependence, but not to a V^-dependence. 

We performed numerical simulations similar to those of Ref. |^ in order to test the analyt- 
ical predictions. The disordered normal region was modeled by a tight-binding Hamiltonian 
on a square lattice (lattice constant a), with a random impurity potential at each site (uni- 
formly distributed between ±^Ud)- The Fermi energy was chosen at Ep = 1.57f/o from the 
band bottom (Uq = fi^ /2ma?'). The length L and width W of the disordered region are 
L = 167a, W = 35a, corresponding to = 15 propagating modes at Ep- The mean free 
path is obtained from the conductance G = {2e'^/h)N{l + L/i)~^ of the normal region in 
the absence of T. The scattering matrix of the normal region was computed numerically at 
e = ±eV, and then substituted into Eq. (H) to obtain the differential conductance. 

In Fig. |1| we show the V^-dependence of G (averaged over some 10^ impurity configu- 
rations), for three values of i. The left panel is for B = 0, the right panel for a flux of 
6h/e through the disordered region. The l^-dependence for i? = is mainly due to the 
contact resistance effect of order N[£/L)'^, and indeed one sees that the amount by which 
G increases depends significantly on The l^-dependence in a T-violating magnetic field 
is entirely due to the weak-localization effect, which should be insensitive to i (as long as 
i/ L <^ 1 <^ N£/L). This is indeed observed in the simulation. Quantitatively, we would 
expect that application of a voltage increases (G) by an amount ^e'^/h for the three curves in 
the right panel, which agrees very well with what is observed. In the absence of a magnetic 
field the analytical calculation predicts a net increase in (G) by 0.79, 0.46, and 0.25 x e^/h 
(from top to bottom), which is again in good agreement with the simulation. 

In conclusion, we have shown that weak localization can coexist with a time-reversal 
symmetry breaking magnetic field in a disordered normal-metal-superconductor junction. 
One needs to apply a magnetic field and to raise the voltage in order to suppress the weak- 
localization correction to the conductance. This leads to an unusual B and \^-dependence 
of the differential conductance, which should be observable in experiments. 

We thank M. J. M. de Jong and J. A. Melsen for help with the numerical simulations. 
This work was supported by the Dutch Science Foundation NWO/FOM. 
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FIGURES 



FIG. 1. Numerical simulation of the voltage dependence of the average differential conductance 
for B = (left panel) and for a flux 6 h/ e through the disordered normal region (right panel) . The 
filled circles arc for an NS junction; the open circles represent the ^-independent conductance in 
the normal state. The three sets of data points correspond, from top to bottom, to i/L = 0.31, 
0.23, and 0.18, respectively. The arrows indicate the theoretically predicted net increase of (G) 
between V = and V » Ec/e. 
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